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Abstract 

We calculate the masses and residues of the triply heavy spin-1/2 baryons using 
the most general form of their interpolating current via QCD sum rules. We compare 
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1 Introduction 



Recently, there have been significant experimental success on the identification and spec- 
troscopy of the heavy baryons containing heavy bottom and charm quarks. By this time, 
all baryons containing a single charm quark predicted by the quark model have been ob- 
served. The heavy A5, S5, S5 and fib baryons with spin-1/2 together with the spin-3/2 

baryon containing a single bottom quark have also been discovered (for a review on the 
spectroscopy of the heavy flavor baryons at the Tevatron see for instance [1]). Recently, the 
CMS Collaboration at CERN reported the observation of the spin-3/2 heavy baryon 
[2]. Following these experimental success on the spectroscopy of the heavy baryons with 
a single heavy quark, the SELEX Collaboration announced their first observation of the 
doubly heavy spin-1/2 S+ baryon with two charm quarks [3-5]. We hope that the LHCb 
detector at CERN will provide us with identification and detection of all doubly heavy and 
triply heavy baryons predicted by the quark model. 

The experimental progresses on the spectroscopy of the heavy baryons have stimulated 
the theoretical works in this respect. There are many works on the spectroscopy of the 
heavy baryons with a single heavy quark in the literature. There are also dozen of works 
dedicated to the spectroscopy of the doubly heavy baryons, however, the number of works 
devoted to the investigation of the properties of the triply heavy baryons are relatively 
few. For some of the theoretical works dedicated to the spectroscopy of the triply heavy 
baryons with different approaches such as the effective field theory, the lattice QCD, the 
QCD bag model, different quark models, the variational approach, the hyper central model, 
the potential non-relativistic QCD and the Regge trajectory ansatz, see for instance [6-18]. 
The masses as well as the masses and residues of the triply heavy baryons for the loffe 
current are calculated within QCD sum rules in [19] and [20], respectively. 

In the present work we extend our previous studies on the spectroscopy and mixing 
angles of the doubly heavy baryons [21-23] to the triply heavy baryons. We calculate the 
masses and residues of the triply heavy spin-1/2 baryons using the most general form of 
the interpolating current via QCD sum rules and compare the results with the predictions 
of the same approach in the case of loffe current [19, 20] as well as predictions of other 
theoretical approaches [6-18]. 

The layout of the article is as follows. In Section 2, we derive QCD sum rules for the 
masses and residues of the triply heavy spin-1/2 baryons. In Section 3, we numerically 
analyze the sum rules for the masses and residues and find the reliable working regions for 
the auxiliary parameters entered to the sum rules. We compare our numerical results with 
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Table 1: The quark contents of the triply heavy spin-1/2 baryons. 

the predictions of the theoretical works existing in the hterature and discuss the obtained 
results. 



2 Masses and residues of the triply heavy spin— 1/2 
baryons 

In order to obtain the QCD sum rules for the masses and residues of the triply heavy 
baryons we start considering the correlation function 

U{q) = i J d'xe"'^{0\r{vQQQ'{x)fjQQQ'{0)}\0), (2.1) 

where tjqqqi is the interpolating current for these baryons and q is their four-momentum. 
The most general form of the interpolating current for the triply heavy spin-1/2 baryons 
can be be written as 

where i, j, k are the color indicies, C is the charge conjugation operator and (5 is an 
arbitrary auxiliary parameter which we shall find its working region. The case /3 = — 1 
corresponds to the loffe current considered in [19, 20]. The heavy Q and Q' quarks for both 
members predicted by the quark model are given in Table 1. From this current one can 
formally obtain the interpolating current of the proton (neutron) by replacing Q ^ u and 
Q' ^ d{Q ^ d and Q' u). 

The correlation function in Eq. (2.1) can be calculated in two different ways. In physical 
or phenomenological side it is calculated in terms of hadronic parts, while in QCD side it is 
evaluated in terms of quarks and gluons. The matching of these two representations then 
gives us the QCD sum rules for physical quantities under consideration. To suppress the 
contributions of the higher states and continuum we apply Borel transformation as well as 
continuum subtraction to both sides of the obtained sum rules. 
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In physical part, by saturating the correlation function with a complete set of hadronic 
states having the same quantum numbers as the interpolating current and by isolating the 
ground state baryons, we get 



— m\ 



(2.3) 



where dots stands for the contribution of the higher states and continuum. The matrix 
element of the interpolating current between the vacuum and the baryonic state is param- 
eterized as 



{^\r]QQQmB{q,s)) ^ \Bu{q,s), 



(2.4) 



where As is the residue of the heavy spin-1/2 baryons and u{q, s) is their Dirac spinor. By 
performing summation over the spins of these baryons, we immediately obtain 

AK^ + ms) 



n(g) = 2"' 

q^ — 



(2.5) 



for the physical side, where we have only two independent Lorentz structures ^ and U to 
calculate the masses and residues. 

In QCD side, the correlation function is calculated using the operator product expansion 
(OPE) in deep Euclidean region. By applying the Wick theorem and contracting out all 
quarks fields, we obtain the following expression in terms of the heavy quarks propagators: 



n(g) 



/ 



Q^'j Qimi 



nkr 



Qlmi 



•n-j Q>'>ni-, nlk , onkr 



0), 

(2.6) 



where S' = CS'^C. 

To proceed in QCD side, we write the coefficients of the selected structures as the 
dispersion integral 



n.(g) = / ^ci., 

s — q^ 



(2.7) 



where Pi{s) are the spectral densities and they are found from the imaginary parts of the 
Ili{q) functions. Here i = 1 and 2 correspond to the structures ^ and U, respectively. Our 
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main task in the following is to calculate these spectral densities. To go further, we need 
the explicit expression of the heavy quark propagator. It is given as 



- Ws 



{2n 



ft + rriQ 
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u 
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where Ki and K2 are the modified Bessel function of the second kind. Using this propagator 
in Eq. (2.6) and after performing lengthy calculations we obtain the spectral densities 



r'tl'max rVmax I 



12(-l + 77)mQmQ.(-l + /3)' 



+ ^2^(3/x^QQ, - 2s) 5 + /3(2 + 5/3) + ^2ml{-l + (3f - UrnqniQ^-l + (3^) 



+ (-1 + ii)r]{ZnQQQ, - 2s) 5 + /3(2 + 5/3) 



+ 



{9IGG) 



2567r4mQmQ/ j^, 
+ 6(-3 + 477) (-1 + + 77)m|,(-l + /3^) + mQUiQ 



d^dr]\ 6(-3 + 4V')(-1 + V' + ^)^q(-1 + Z^^) 

48V'^(l + ^^)+V' -63 



+ 68?7 - 30/3 + 877^9 + (-63 + ^%r])(3'^ + 2(-l + r/) - 3 3 + /3(2 + 3/3) 



+ 2r/ 



5 + /3(2 + 5/3) 



(2.9) 
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and 



*^ V'mm *^ Vrrtirt I 



r){-l + + ri)mQ,{nQQQi - s)(-l + (5f 



+ 6^(-l + + v)mQ{fXQQQ' - s){-l + P^) + mlrriQ, 5 + /3(2 + 5/3) 

{glGG) /•V'ma:. /-r^ma:, ^^^^ 
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-2(-l + 77)77mQm^,(-l + ^)^ 



- 2^3^(-l + /?)[- 9mQ.(/iQQQ. - + /3) + r]{2iJiQQQ, - 3s) (mQ(-l + /3) 



+ QmQ,{l + /5))] + ^PniQ (^rf{iiQQQ, - s){-l + + 2mQmQ,(-l + 



+ 77 -5m^,(-l + ^)2-2mQmQ,(-l + ^')+m| 5 + ^(2 + 5/3) 
+ V'' - 4mJmQ,(-l + ^2) + 77'(7/iQQQ' - 9s)(-l + ^) [mQ(-l + ^) + 6mQ,(l + ^)] 
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2ri%2f,QQQ, - 3s)(-l + /3) ruQi-l + /5) + 6mQ.(l + /3) 



-1 + /32) + 12mQm|,(l + Z^^) - m| 5 + /3(2 + 5/3) 



18mQ/(/iQQQ 
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(2.10) 



where, 
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I, - Ami - ^(s + - 4m|)2 - 4mJ,s 
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(2.11) 



Now we match the two sides of the correlation function for each structure to get the 
QCD sum rules for the masses and residues of the triply heavy baryons. After applying 



the Borel transformation and continuum subtraction to suppress the contributions coming 
from the higher states and continuum, we get 



7(2mQ+mQ,)2 

-ml r^o 
i2 — e. 



X^niBe — I dsp2{s)eM^ , (2-12) 

J(2mQ+mg,)2 

where and Sq ^^I'c Borel mass parameter and continuum threshold, respectively. By 
eliminating the residues from the above equations, we can calculate the masses from either 

m% = z = 1 or 2, (2-13) 



or 



— s 



!lL^^m^,Ydsp,{s)ei 
ruB = — —■ (2-14) 

3 Numerical results 

Now we are ready to numerically analyze the obtained sum rules in the previous sec- 
tion and calculate the numerical values of the masses and residues of the triply spin- 
1/2 heavy baryons. For this aim we take the quark masses as both their pole values 



mi, = (4.8 ± 0.1) GeV and m^ = (1.46 ± 0.05) GeV [24] as well as their MS values 
ffibimb) = (4.16 ±0.03) GeV, rhcirhc) = (1.28 ±0.03) GeV [25]. For the numerical value of 
the gluon condensate we use (g^.GG) = 47r2(0.012 ± 0.004) GeV^ [24]. 

The sum rules obtained in the previous section incorporate also three auxiliary param- 
eters which we find their working regions in the following. These parameters are the Borel 
mass parameter M^, the continuum threshold Sq and the general parameter (3 enrolled to 
the general current of the baryons under consideration. The working regions of these pa- 
rameters are found such that the variations of the masses and residues are very weak with 
respect to them. 

The continuum threshold Sq is not completely capricious but it is related to the energy 
of the first excited state. We have not adequate knowledge about the first excited states of 
the baryons under consideration, but our analysis show that when we choose the continuum 
threshold in the intervals Sq = (140 - 148) GeV'^ and Sq = (74 - 81) GeV'^ respectively for 
the ftbbc and flccb baryons in the case of pole quark masses, the results very weakly depend 
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on So- In the case of MS values of the quark masses, the working regions for the continuum 
threshold are obtained as Sq — (117 — 125) GeV"^ and Sq — (64 — 70) GeV"^ for the baryons 
flbbc and ^ccb, respectively. 

Now we proceed to find the working region for the Borel mass parameter M^. The upper 
bound on this parameter is found demanding that the pole contribution is high compared 
to the contributions of the continuum and higher states. This means that the relation, 

) 1/2, (3.15) 

J {ImQ+rriQiY 

is satisfied. This condition leads to the following upper values for M^: 



J (2m, 



{2mQ+mQ 



, 22 GeV^, for Q^bc 
MLx-< ' 3.16 

^ 18 GeV^, for Q.^- 

The lower bound on is calculated requiring that the contribution of the perturbative 
part exceeds the nonperturbative contributions. From this restriction we obtain 

„ f 12 GeV^, for Qbbc 
MLn ={ . 3.17 

\ 9 GeV^ for Q,,b- 

Our final task is to determine the working region for the auxiliary parameter (3. Instead 
of discussing the variations of the physical observables with respect to this parameter in 
the interval (— oo,+oo), we define (3 — tanO and look for the variations with respect to 
cos9 in the interval [—1, 1]. Our numerical results show that in the domains [—0.5, —0.9] 
and [0.5,0.9], the residues depend weakly on cosO. Here we should mention that the loffe 
current corresponds to cosO — —0.71 and lies inside the rehable region. Note also that as 
the mass sum rule is the ratio of two expressions including ^ in Eq. (2.13), the masses show 
a very good stability with respect to cos9 in the whole allowed region. 

Considering the working regions of the auxiliary parameters we obtain the numerical 
values for the masses and residues of the triply heavy spin-1/2 baryons as presented in 
Tables 2 and 3 for both structures. For comparison we also present the numerical predictions 
of other theoretical approaches such as the modified bag model [12], the relativistic three- 
quark model [13], the non- relativistic three-quark model [14] and QCD sum rules but loffe 
current [19,20]. As far as the masses are considered, our central value results are slightly 
higher than the other predictions. The closest results to our predictions are the results of 
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This work (^) 


This work (U) 


[20] 


[12] 


[13] 


[19] 


[14] 


^bbc 


11.73 ±0.16 


11.71 ±0.16 


11.50 ±0.11 


11.139 


11.280 


10.30 ±0.10 


11.535 


^ccb 


8.50 ±0.12 


8.48 ±0.12 


8.23 ±0.13 


7.984 


8.018 


7.41 ± 0.13 


8.245 


^bbc 


10.59 ±0.14 


10.56 ±0.14 


10.47 ±0.12 










^ccb 


7.79 ±0.11 


7.74 ±0.11 


7.61 ±0.13 











Table 2: The masses of the triply heavy spin-1/2 baryons in the units of GeV compared 
with other theoretical predictions. For the baryons with over-line, the MS values of the 
quark masses have been used. 





This work (^) 


This work (U) 


[20] 


^bbc 


0.53 ±0.17 


0.45 ±0.15 


0.68 ±0.15 


^ccb 


0.38 ±0.13 


0.30 ±0.10 


0.47 ±0.10 


^bbc 


0.85 ±0.28 


0.65 ±0.22 


0.68 ±0.15 


^ccb 


0.56 ±0.18 


0.38 ±0.13 


0.47 ±0.10 



Table 3: The residues of the triply heavy spin-1/2 baryons in the units of GeV^ compared 
with the existing theoretical prediction for loffe current. For the baryons with over-line, 
the MS values of the quark masses have been used. 

the non-relativistic three-quark model [14] and QCD sum rules with the loffe current [20], 
respectively. The lower predictions for the masses belong, respectively, to QCD sum rules 
with the loffe current [19] and the modified bag model [12]. From Table 2 we see that the 
two structures in our case give approximately the same results. This Table also shows that 
the results depend on the quarks masses considerably and change (9-10)% when going from 
the pole masses to the MS scheme masses. Here, we should mention that considering Eq. 
(2.14) does not affect considerably the results of masses presented in Table 2. 

In the case of the residues, against the [20] our results depend on the quark masses 
such that when going from the pole to the MS the results change (21-37)%. This is an 
expected results that the residues depend more on quark masses in comparison with the 
baryon masses. From Table 3 it is also clear that the results depend on the choice of the 
structure. The structure U give the results (15-30)% lower than those of the structure ^. In 
the case of the pole masses of the quarks, our predictions on the residues are considerably 
small in comparison with those of the [20]. The maximum difference between two works 
belongs to the Jlccb baryon and U structure which is approximately 36%. For the residues 
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in MS scheme, our result is very close to that of the [20] for the structure U and flbbc, while 
the maximum difference with 23% between predictions of two studies belongs to the flccb 
baryon and also the structure U. 

In conclusion, we have calculated the masses and residues of the triply heavy spin-1/2 
baryons using the most general form of their interpolating current in the framework of 
QCD sum rules. We found the reliable working regions for each of the auxiliary parameters 
entered the calculations. Our results on the masses are slightly higher than the previous 
predictions like the modified bag model, different relativistic and non-relativistic quark 
models as well as QCD sum rules for the loffe current. We obtained residues considerably 
different than the existing predictions via QCD sum rules with loffe current. We hope that 
the LHC at CERN will provide opportunity to experimental studying of these baryons in 
near future. 
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